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李 建 春 , 女 ,博士 ,东南 大 学 首席 教授 .博士 生 导 师 , 国 家 杰出 青年 科 

学 基金 获得 者 .江苏 省 双 创 人 才 以 及 江苏 省 双 创 团队 负责 人 。 现 任 《 Geo- 

mechanics and Geophysics for Geo — Energy and Geo - Resource) 联合 主编 、 

《Rock Mechanics and Rock Engineering》 副 主编 ,同时 担任 国际 岩石 力学 与 
岩石 工程 学 会 (ISRM) 岩石 动力 学 专 委 会 主席 。 

长 期 致力 于 岩石 动力 学 应 用 基础 研究 ,开展 了 应 力 波 通过 岩石 节理 

的 连续 -- 非 连续 理论 方法 研究 , 提出 了 应 力 波 通过 节理 的 二 维 时 域 递 归 

分 析 方法 ,建立 了 充填 节理 的 薄弱 层 界 面 模 型 和 岩 体 夭 弹 性 动态 等 效 连 续 介质 模型 ,揭示 了 应 

. 力 波 在 岩 体 中 传播 规律 。 近 年 来 ,作为 负责 人 主持 国家 杰出 青年 科学 基金 项 目 、 国 家 自然 科学 

P 基金 重点 项 目 和 重点 国际 合作 项 目 、 国 家 重点 研发 专项 课题 。 在 国内 外 权威 期 刊 发 表 论文 100 


) 余 篇 ,其 中 SCI 收录 论文 80 余 篇 ,出 版 英文 专著 2 部 中文 专著 1 部 . 参 编 英文 专著 2 部 。 
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So 岩 体 中 应 力 波 传播 规律 研究 方法 进展 
pt 
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OF 要 :应 力 波 在 岩 体 中 传播 规律 的 研究 是 分 析 工 程 岩 体 动力 稳定 性 的 基础 。 天 然 岩 体 中 包含 的 
大 量 节理 对 应 力 波 的 传播 规律 具有 显著 的 影响 。 应 力 波 在 岩 体 中 传播 时 会 出 现 幅 值 的 衰减 以 及 波 
加 的 降低 ,在 此 过 程 中 ,节理 在 应 力 波 的 作用 下 会 出 现 张 开 、 闭 合 和 滑 移 变形 。 目 前 ,应 力 波 在 岩 体 
中 的 一 维和 二 维 传播 规律 分 析 方法 比较 成 熟 ,这 些 方法 主要 分 为 非 连续 介质 方法 与 动态 等 效 连 续 
介质 方法 以 及 双 尺 度 非 连续 褐 合 方法 。 借 助 这 些 方法 ,研究 者 们 开展 了 线性 和 非 线性 节理 对 应 力 
波 传播 规律 影响 的 研究 。 本 综述 将 详细 地 总 结 上 述 3 种 方法 ,探讨 它们 的 优 缺点 ,并 提出 值得 进 一 
步 研 究 的 问题 。 
关键 词 : 岩 体 ;应 力 波 ;节理 ;位 移 不 连续 方法 ;等 效 介质 方法 
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Progress of methods for wave propagation across rock masses 


LI Jianchun' , FAN Lifeng’ , LI Zhengliang! 
(1. Engineering Research Center of Safety and Protection of Explosion & Impact of Ministry of Education, 


Institute of Future Underground Space , Southeast university ,211189 Nanjing, China; 


2. College of Architecture and Civil Engineering , Beijing University of Technology , 100124 Beijing, China) 


Abstract ; Studying stress wave propagation across jointed rock masses is essential for the stability of rock 
engineering. Rock masses in nature include many joints , which have great influence on stress wave propa- 
gation. During wave propagation across rock masses, both of the amplitude and speed of the wave attenu- 
ate. Meanwhile, the interaction between the stress wave and the joints usually results the joints open, close 
and slippery. Currently, some analytical methods have been successfully applied-for studying one-dimen- 
sional or two-dimensional wave propagation across rock joints. These analytical methods are divided into 


three types: the discontinuous medium method, the equivalent continuous medium methods and double- 


^ scale discontinuous method. Based on these two methods, the effect of linear and nonlinear joints on wave 


“propagation has been studied effectively. The present paper is to introduce the discontinuous and equiva- 


lent continuous media methods, discuss their applicability , and suggest the further problems to study. 


Key words:rock mass; stress wave; joints; displacement discontinuity method; equivalent continuous 
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j 岩 体 中 应 力 波 传播 规律 的 研究 是 目前 国内 外 研 


完 的 热点 问题 之 一 ,也 是 岩石 力学 的 难点 问题 ,该 研 
窒 绪 果 对 分 析 和 评估 动态 扰动 下 岩 土工 程 的 安全 稳 
定 有 重要 意义 。 例 如 ,爆破 引起 的 振动 通常 以 应 力 
波 的 形式 在 岩 体 中 传播 ,并 对 周围 环境 造成 影响 。 
对 振动 的 预测 与 控制 是 确保 爆破 工程 顺利 开展 的 基 
础 > 研究 应 力 波 在 岩 体 中 传播 规律 可 以 为 爆破 振动 
的 洪 确 预测 提供 理论 依据 。 因 此 深入 理解 和 研究 应 
力 波 在 岩 体 中 的 传播 规律 意义 重大 。 

“由 于 岩 体 中 大 量 节理 的 存在 ,应 力 波 在 传播 过 
程 中 频率 和 幅 值 都 会 发 生 衰减 。 应 力 波 在 节理 处 的 
反射 .折射 和 散射 是 其 衰减 的 重要 因素 "” 。 不 同 尺 
度 的 结构 面 对 应 力 波 传播 规律 具有 不 同 的 作用 机 
MD 。 为 此 ,国内 外 学 者 主要 针对 岩石 中 的 裂 颖 和 
岩 体 中 的 节理 开展 了 理论 分 析 , 得 到 应 力 波 通过 着 
石 或 岩 体 时 发 生 的 衰减 和 频 散 ”” 。 岩 体 中 包含 很 
多 类 型 的 不 连续 结构 面 , 按 其 尺度 可 分 为 五 级 ,其 中 
IL AN IV 级 结构 面 为 呈现 不 同 延 伸 性 的 节理 ,在 
岩 质 隧道 . 边 坡 和 基础 工程 中 广泛 存在 ,影响 着 岩 体 
结构 完整 性 和 物理 力学 性 质 * 。 从 工程 角度 来 说 ， 
节理 对 应 力 波 传播 规律 的 影响 不 容 忽视 。 现 有 的 理 
论 分 析 方 法 主要 分 为 两 大 类 , 即 非 连 续 介 质 方法 与 
动态 等 效 连续 介质 方法 。 这 些 方法 都 有 各 自 的 优点 
和 局 限 性 。 因 此 ,有 必要 对 现 有 的 方法 进行 详细 地 


投稿 网 站 : http ;// cjam. xjtu. edu. cn. 微 信 公众 号 :应 


综述 ,发 现 它 们 的 适用 范围 ,从 而 为 将 来 进一步 的 研 
提供 指导 。 


1 非 连续 介质 分 析 方 法 


非 连续 介质 方法 是 分 析 应 力 波 在 岩 体 中 传播 规 
律 的 有 效 方法 之 一 。 其 代表 方法 是 位 移 不 连续 方 
法 ,该 方法 将 节理 视 为 黏 结 ( 节理 面 的 位 移 连 续 且 作 
用 在 节理 前 后 的 应 力 连续 ) AN SEE ht CBE TEE 
位 移 不 连续 但 作用 在 节理 前 后 的 应 力 连续 ) 或 完全 
张 开 ( 位 移 不 连续 但 作用 在 节理 上 应 力 为 零 ) 的 边 
界 条 件 。 该 方法 可 用 于 分 析 应 力 波 垂直 或 者 以 任意 
角度 斜 人 射 单条 或 一 组 节理 。 


1.1 应 力 波 通过 单条 节理 的 传播 


亚 级 岩 体 中 包含 着 大 量 的 区 域 性 节理 ,这 类 结 
构 面 的 特征 尺寸 相 较 于 波长 具有 延伸 好 但 厚度 薄 的 
特点 。 在 假设 节理 为 无 限 延伸 的 情况 下 , 现 有 的 研 
究 对 应 力 波 通过 单条 节理 开展 了 大 量 的 研究 。 
Schoenberg“ 将 节理 视 为 应 力 连 续 、 位 移 不 连续 的 
线性 滑 移 界面 ,计算 了 简 谐 波 以 任意 角度 入 射 节理 
时 的 透射 和 反射 系数 。 当 任意 压缩 波 (P UE) AMETE 
介质 1 和 2 之 间 的 界面 时 ,应 力 波 会 发 生 透 射 和 反 
射 (图 1)。 
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介质 2 


图 1 平面 简 谐 P 波 在 介 质 1 和 2 之 间 
非 连续 界面 上 的 透射 和 反射 


Fig.1 Reflection and transmission of a plane harmonic P 


wave on the interface between media 1 and 2 


[ —k,cos6,, [sing0， 一 上 .cosg + ez, C0826. 

à E 

m E A 

Kk, sing, -k,cos0,  k,sinĝ y — iw sin20,, 

SN Zp2 

4C0820, — z4sin20, 2,5 C0820,; 
2 
: £o . 

C sin20,, z 0820, - sin20,, 

[人 


其 中 :R,R.、T, 和 7 表示 反射 系数 和 透射 系数 , 即 
hee 
Ego W P OIH S W) ES 质 1 中 的 波 阻 抗 ， 
Zo? = Prey 301 为 介质 1 的 密度 ie Cus AD PE BF 
ASE) 在 介质 1 中 的 波 速 ;P 波 和 S 波 的 传播 广 
m E a 


sin, _ Sin, _ sind _ sing P 
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图 2 节理 法 向 力学 特性 的 双 曲 线 模型 示意 网 


Fig.2  Hyperbolic elastic model of normal deformation 


behavior of a joint 
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如 果 该 界面 为 节理 ,具有 一 定 的 法 向 和 切 向 刚 
度 ,根据 位 移 不 连续 假设 ,作用 在 节理 上 的 应 力 连 续 
但 是 节理 前 后 的 位 移 不 连续 , 即 
04 = 0 = g. 
Tia 一 Tu = Tus 1 
Wy —- us, =0,/k, i 
Uy 一 La ETa k, 
HP :0, 为 作用 在 节理 上 的 法 向 应 力 ; 下 标 1 和 2 对 应 
介质 1 和 2;7,. 为 作用 在 节理 上 的 切 疝 应 力 ;u, 和 ww, R 
示 和 z [a] LAB Sk, 和 不 为 节理 的 切 向 和 法 向 刚度 。 
将 入 射 波 、 反射 波 和 透射 波 引起 的 位 移 、 应 力 代 入 式 
(1) 可 得 到 关于 透射 和 反射 系数 的 矩阵 方程 ”为 


k. sin, — iwz,sin26., 


-R [| 一 k.cosQ,, 
. p ; 
k.cos0,, — iwz..c0s26., R k,sin6,, 
s 
2 x r |7 z,, 0820. (2) 
7 4981 Ao p 2 
z 
sb. 
LT, - sin20,, 
Z C0S20 2 L Zp 


节理 中 充填 的 水 会 引起 节理 的 黏 性 性 质 。 因 此 ,他 
们 推导 了 频 域 内 考虑 节理 黏 弹性 特性 的 应 力 波 反 
射 .透射 系数 的 理论 公式 ,并 通过 试验 揭示 了 超声 波 
通过 完整 饱和 岩石 试 样 时 透射 S 波 脉 冲 与 超声 波 通 
过 包含 节理 的 岩石 试 样 时 透射 $ 波 脉 冲 之 间 的 差 
异 。 除 此 之 外 ,节理 中 的 填充 物 也 可 能 导致 节理 具 
i BREE?! ,Zhu SE EAA 
BUREAU FE, EAB E.A TY LY Bh 
弹性 特性 对 透射 波 频谱 特性 的 影响 。 

然而 ,在 处 理 爆炸 所 产生 的 应 力 波 传播 问题 时 ， 
研究 人 员 发 现 当 应 力 波 的 幅 值 较 高 时 应 力 波 会 引起 
节理 的 非 线 性 变形 行为 2 7, Miller 4:7? 推导 了 平 
面 简 谐 剪 切 波 垂直 通过 具有 剪 切 滑 移 特性 节理 的 透 
射 和 反射 系数 ,此 时 ,节理 面 上 的 剪 切 应 力 为 滑 移 量 
和 滑 移 速率 的 非 线性 函数 。 由 于 难以 获得 非 线 性 方 
程 的 解析 解 , Miller 等 采用 了 将 非 线 性 问题 简化 为 线 
性 问题 的 近似 分 析 方法 开展 研究 ,结果 表明 透射 系 
数 取决 于 入 射 应 力 波 的 振幅 与 节理 的 抗 剪 强度 。 

相 比 于 频 域 下 的 分 析 方 法 ,时 域 下 时 间 差 分 分 
析 方 法 能 更 好 地 适用 于 节理 非 线 性 响应 的 研究 。 
Zhao 等 "采用 改进 的 特征 线 法 分 析 了 平面 P YE 


Schoenberg! 5! 发 现 即使 假设 节理 刚度 与 频率 无 
关 , 反 射 和 透射 系数 仍 依赖 于 人 射 波 的 频率 。 卢 文 
波 中 通过 开展 应 力 波 与 线性 可 滑 移 岩 石 节理 之 间 相互 
作用 的 研究 ,揭示 了 节理 对 应 力 波 的 高 频 滤波 效应 。 

当 岩 体 处 于 饱和 状态 时 ,Pyrak-nolte 等 ' 中 认为 


直通 过 单条 节理 的 透 反 射 情况 ,其 中 ,假定 节理 具有 
Barton 等 ! E ] 提 出 的 双 曲 线 型 (B-B 模型 ) 法 向 非 线 
性 行为 (图 2)。 该 B-B 模型 表示 为 

kad, 


RTA (4) 
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其 中 :d, 和 di 是 广 理 的 法 向 闭合 量 和 最 大 允许 闭 
合 量 ;为 节理 的 法 向 初始 刚度 。Barton 等 ” 指出 
d 和 各 可 根据 实验 室 测试 得 到 的 节理 面 粗糙 度 
(JRC) 节理 面 抗 压强 度 (JCS) 和 节理 的 平均 缝 宽 
(ai) 来 确定 。 节 理 处 应 力 波 的 透 反 射 及 其 特征 线 如 
图 3 所 示 。 


节理 
l 
入 射 波 一 | 


三 一 透射 波 


反射 波 < 


振 速 边界 条 件 


(a) 节理 处 应 力 波 的 透 反 射 


节理 


v) | veo 
eG) | eoim 


pae o x, o(0)0 b x 
(b) 节理 处 应 力 波 的 特征 线 
图 3 节理 处 应 力 波 透 反 射 及 其 特征 线 


g.3 Stress wave normally transmission across a joint, 


and the corresponding left-and right-running 


= : characteristics in the x-t plane 

二 特征 线 方法 认为 节理 为 非 连续 边界 ,此 时 ,作用 
在 节理 上 的 应 力 和 节理 面前 后 的 位 移 应 满足 如 下 
RPI, 


o (x,,t) =o (x,t) (5) 
Go nu Gon SETTE] 
考虑 节理 面前 后 质点 振动 速度 , 式 (6) 可 进 一 
步 表 示 为 
v (x,,t) -v'(x,,t) 
| 90 (x, ,t) i ka (7) 
at Us + Lon 07d] F 
特征 线 方法 认为 沿 右 行 和 左 行 特征 线 以 及 应 力 
与 质点 振动 速度 的 关系 为 


zo(x,1) +o(x,t) = const, AFTER (8) 
20(x,t) a (x,0) = const, 左 行 特征 线 (9) 


其 中 z 为 应 力 波 的 阻抗 。 当 初始 时 刻 节 理 面前 后 质 
点 振动 速度 、 节 理应 力 以 及 边界 处 的 入 射 波 已 知 时 ， 
由 特征 线 方法 可 得 到 节理 面前 后 质点 振动 速度 和 节 
理应 力 。 应 该 注意 的 是 , 当 采 用 特征 线 求解 时 , 式 
(7) 转 化 为 时 间 差 分 方程 。 

国内 外 科研 工作 者 基于 特征 线 方 法 开展 了 大 量 
的 应 力 波 垂 直入 射 单条 或 一 组 平行 节理 的 传播 规律 
研究 。 例 如 ,Zhao 等 计算 了 剪 切 波 垂直 人 射 
具有 库伦 滑 移 特性 节理 的 透射 系数 (图 4) ,结果 表 
明 入 射 波 的 前 应 力 峰值 与 节理 的 剪 切 强度 之 比 是 影 
响应 力 波 透射 特征 的 最 重要 因素 。 在 此 基础 上 , 俞 
络 等 中 提出 了 节理 的 三 参数 非 线 性 本 构 模 型 ,并 
基于 该 模型 采用 特征 线 法 探讨 了 了 波 通过 单条 非 线 
性 节理 的 透射 情况 。Fan °° 分析 了 具有 不 同 加 钾 
载 变形 特性 的 非 线性 节理 对 应 力 波 传播 规律 的 影 
Wi], Li 等 2 基于 特征 线 方 法 分 析 了 了 波 垂直 入 射 
时 节理 的 张 开 闭 合力 学 特性 对 透射 波 波形 、 幅 值 和 
频谱 的 影响 ,并 通过 分 离 式 霍 普 金森 (SHPB ) 冲击 实 
验 验 证 该 理论 分 析 结 果 的 正确 性 。Fan 45 EF g- 
A 非 线性 节理 本 构 模 型 研究 了 了 波 通过 节理 贿 体 时 
的 幅 值 衰减 和 速度 降低 效应 。Fan 等 分 析 了 软 硬 
岩 之 间 非 线性 结构 面 对 应 力 波 传播 规律 的 影响 。 特 
征 线 方 法 能 够 快速 便捷 地 分 析 应 力 波 通过 节理 岩石 
的 一 维 波 传播 规律 问题 ,但 不 适用 于 应 力 波 和 斜 人 射 
节理 的 情况 。 


1 
1 
1 
1 
1 
1 
1 
1 
1 
P 
ik, | 
Pa 
$00. 
1 
1 
1 
1 
d, 


图 4 节理 的 切 向 库伦 滑 移 模型 

Fig.4 Coulomb slip model of a joint 

频 域 分 析 方 法 ,如 公式 (1) ,可 用 于 分 析 应 力 波 

以 任意 角度 ,特别 是 斜 入 射 节 理 的 情况 。 由 于 频 域 

分 析 方法 只 针对 岩 体 为 线性 系统 情况 , 即 节理 和 其 

周围 的 岩石 都 是 线 弹性 ,并 不 适用 于 节理 非 线性 力 

学 行为 的 情况 ,因此 对 于 非 线 性 节理 ,就 必须 寻求 更 
合适 的 方法 来 研究 应 力 波 斜 入 射 节 理 问 题 。 

当 应 力 波 以 任意 角度 斜 人 射 节理 时 ,为 了 分 析 

节理 非 线性 对 应 力 波 传播 规律 的 影响 ,Li 等 ”提出 
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了 应 力 波 二 维 传播 规律 的 时 域 递 归 分 析 方 法 。 在 该 
ro ce 
成 的 单元 体 满足 力 的 平衡 条 件 , 且 波 阵 面 上 满足 
basins glee ed spe 
情况 ,时 域 递 归 分 析 方法 从 分 析 应 力 波 和 节理 之 间 
的 相互 作用 出 发 (如 图 5 所 示 ) ,通过 分 别 分 析 节 理 
前 后 单元 体 应 力 状 态 ,得 到 节理 面 的 质点 振动 速度 
以 及 作用 在 节理 上 的 应 力 ™1。 
0 =2,c0s26,(v,, Vrp) — z,sin20,vg, 


z,sin26,tand, 


T = t 0 ( Vip z Un, ) ~ z,C0820, Vps 
an, 
" (10) 
0 =2,c0s26.v,, +Z sin20,vy, 
+ _2,sin26,tand, 
T mW, — 2,00820,v4, 
tang E Se 
p 
o; = 0086, (iy -om) + sin ny, 
v. =sing, ( Uy + Upp) + cosh Vg, (11) 


v. — cOSÜ Vr + Sinh Vy 


€ v; =sing pVTp T COSO Vrs 

二 如果 节 理 的 法 向 力学 特性 表现 为 B-B 模型 (图 
"A DIR] AERE EOS Coulomb WEBEBERI CRI 4) , Du 
u, —u, -g/(k, +o/d,,,,) 
u` -u =7/k., d, xd, (12) 
kir T= +T, = tglang, d, » d, 


TU k, 表示 节理 的 剪 切 刚 度 ;r, 和 d, 分 别 为 节理 


节理 
(d) 反射 S 波 
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(e) 透射 P 波 


AdE HHI 
iv CI 1 FAT] | I) 849 


的 剪 切 强 度 和 对 应 的 剪 切 变 形 ;p 为 节理 的 摩擦 角 。 
式 (12) 可 改 成 质点 振 速 差分 形式 ,考虑 节理 面 的 位 
移 不 连续 而 应 力 连续 条 件 , 将 式 (10) ~ (11) FLASK 
(12) ,可 得 到 应 力 波 以 任意 角度 斜 人 射 单条 节理 的 
传播 方程 为 


Vp, (t£) 4d 4 va (t) 

--B Av B C 13 
bo ind el 
[rE] eap, p tgo D] 
vg (t+ At) i Up. (t) 

1 vt) 
G F 14 
b MM 


其 中 :4 ~F 为 参数 矩阵 ,具体 表达 式 见 文献 [40 ] 。 
当 系 统 的 初始 条 件 和 边界 条 件 已 知 时 ,基于 式 (13) 
开展 计算 就 可 得 到 节理 处 随时 间 变 化 的 透射 波 和 反 
射 波 。Li 等 ("| 采用 时 域 递归 分 析 方 法 研究 了 任 
意 波形 应 力 波 通过 线性 和 非 线性 节理 的 传播 规律 。 
Zou 等 “基于 时 域 递归 分 析 方 法 开展 了 从 数值 模 
拟 或 平面 波动 平台 获得 的 数据 中 分 离 节 理 面 前 后 
透 . 反 射 P 波 ,SV 波 的 方法 。Li EO 基于 时 域 递归 
分 析 方 法 分 析 了 节理 的 张 开 闭 合力 学 特性 对 了 波 斜 
入 射 节理 后 能 量 衰减 的 影响 。Liu 等 “1 利用 时 域 递 
归 分 析 方 法 分 析 了 不 同 应 力 下 应 力 波 以 任意 角度 斜 
入 射 非 线 性 节理 的 传播 规律 。 汪 书 敏 等 ”| 基于 时 
域 递归 分 析 方 法 分 析 了 节理 黏 弹性 特性 对 应 力 波 传 
播 规 律 的 影响 。 


节理 节理 
(b) 入 射 P 波 


(c) 反射 P 波 


节理 节理 


(D 透射 S 波 


图 5 时 域 分 析 方 法 的 力学 模型 (针对 了 波 斜 人 射 情况 ) 


Fig.5 Stress analysis for an P wave incident on a joint 
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1.2 应 力 波 通过 一 组 平行 节理 的 传播 


自然 界 的 节理 通常 成 组 出 现 ,节理 在 每 个 节理 
组 中 基本 呈现 为 保持 一 定 距 离 的 平行 状态 (1 。 当 
相 邻 节理 之 间距 离 较 大 时 ,透射 波 主要 取决 于 首 达 
波 ,从 而 可 忽略 应 力 波 在 节理 之 间 的 多 次 透 反射 产 
生 的 次 达 波 。Pyrak-Nolte 等 '” 基 于 忽略 节理 间 多 
次 透 反射 的 假设 ,提出 了 一 个 计算 应 力 波 垂直 通过 
多 条 平行 节理 后 透射 系数 的 简化 公式 : [Ty | = 
|, I HEHE |7, | 为 应 力 波 垂直 透 过 单条 节理 时 透 
射 波 与 人 射 波 的 幅 值 之 比 ; N 为 节理 数量 。 该 简化 
公式 认为 随 着 应 力 波 通过 节理 数量 的 增加 ,应 力 波 
透射 系数 成 寡 次 方 下 降 。 然 而 当 节理 的 间距 较 小 时 
( 相 较 于 波长 ) ,多 次 透 反 射 的 影响 增强 ,该 公式 将 
不 躲 适用 。 为 了 考虑 应 力 波 在 节理 间 多 次 透 反射 的 
骂 油 ,近年 来 一 些 新 的 理论 分 析 方法 也 不 断 涌 现 ,如 
dtc sea RAE MA, 

gi He f 3 D UE E T s b OR APE AL BL, Zhao 
A pisse et T DIOE P ae SV 流 
以 在 意 角度 斜 入射 多 条 平行 的 线性 节理 的 衰减 规 
向 ; 间 开 展 了 影响 参数 分 析 , 如 节理 的 间距 、 数 量 和 
网 岩 以 及 入 射 角 等 。Perino 等 "使 用 矩阵 传递 法 分 
BEC Ta SH 波 倾斜 人 射 非 线性 节理 的 情况 ,并 采 
Ji: Coulomb 滑 移 模型 描述 节理 的 切 向 变形 行为 。 
Huang °°) 用 矩阵 传递 法 计算 了 应 力 波 通过 具有 
条 弹性 的 一 组 平行 节理 的 透射 和 反射 系数 。 值 得 注 
意 鲍 是 ,矩阵 传递 法 获得 了 频 域 解 ,但 是 对 于 节理 为 
韭 线性 的 问题 , 它 需 要 将 非 线 性 问题 简化 为 线性 问 
题 洲 能 得 到 非 线 性 方程 的 近似 解 。 

除了 矩 阵 传递 法 , Cai 等 由 改进 了 特征 线 法 来 
考虑 应 力 波 在 节理 之 间 的 多 次 透 反 射 对 透射 波 的 波 
形 和 能 量 的 影响 。Zhao 403 利用 改进 的 特征 线 
法 计算 了 P 波 .SV 波 垂直 通过 一 组 线性 或 非 线性 平 
行 节理 后 的 透射 系数 。 在 采用 特征 线 法 来 分 析 问 题 
时 ,需要 首先 将 微分 方程 转化 为 时 间 差 分 方程 ,由 此 
得 到 用 质点 振动 速度 和 应 力 表达 的 应 力 波 传播 方 
程 ,然后 开展 数值 求解 。 特 征 线 方法 可 以 对 具有 任 
意 波形 的 人 射 波 进行 直接 求解 分 析 。 需 要 注意 的 
是 ,特征 线 法 要 求 计算 的 时 间 步 足够 小 ,以 保证 计算 
的 准确 性 ,此 外 ,特征 线 法 难以 处 理应 力 波 以 任意 角 
度 斜 人 射 的 情况 。 

Li 等 55 提 出 的 虚拟 波源 法 也 被 用 于 考虑 应 力 

波 在 节理 之 间 的 多 次 透 反 射 对 透射 波 的 影响 。 该 方 
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法 中 ,虚拟 波源 被 放置 于 节理 处 , 当 任 何 透射 或 者 反 
射 波 到 达 该 节理 位 置 时 虚拟 波源 都 会 产生 新 的 应 力 
波 。 因 此 ,应 力 波 在 岩 体 中 传播 一 段 距离 后 所 产生 
的 透射 波 可 以 看 作 是 应 力 波 通过 该 段 节理 的 直接 传 
播 ( 首 达 波 ) 及 由 虚拟 波源 产生 的 透射 波 ( 次 达 波 ) 
HII, Zhu 等 5 将 虚拟 波源 的 概念 引入 到 位 移 
不 连续 方法 ,分 析 了 应 力 波 垂直 人 射 一 组 线性 节理 
的 情况 。 相 比 于 特征 线 法 ,Zhu 5&0" 认为 虚拟 波源 
法 具有 以 下 优点 :首先 ,虚拟 波源 法 可 以 获得 准确 的 
理论 解 ,而 特征 线 法 需要 将 微分 方程 转化 为 差分 方 
程 ,由 此 不 可 避免 地 产生 计算 误差 ;其 次 ,虚拟 波源 
法 可 以 轻松 分 离 不 同时 刻 到 达 的 透射 波 ; 此 外 ,虚拟 
波源 法 可 以 扩展 到 研究 更 复杂 情况 的 波 传 播 问 题 ， 
例如 应 力 波 斜 信 射 节理 的 情况 "。 虚 拟 波源 法 的 
优点 非常 明显 ,但 是 也 存在 它 的 局 限 性 , 即 难以 处 理 
节理 非 线 性 的 问题 。 

当 应 力 波 在 相 邻 节理 之 间 岩 石 中 的 传播 规律 得 
到 后 ,时 域 递 归 分 析 方 法 也 可 用 于 分 析 应 力 波 斜 人 
射 一 组 线性 和 非 线 性 平行 节理 的 情况 i。 如 果 
相 邻 节理 之 间 岩 石 为 弹性 介质 , 则 岩石 中 的 应 力 波 
为 时 间 平 移 函 数 , 即 

) a5 


D. [77]0- S 

vy, (4) Uy, jel cosÓ, '6 

v. (t) DR 

| i -| : lia | (16) 
vy (t) cos, * c, 


其 中 :5 为 节理 间距 ;J 为 节理 编号 ;下 标 Ip .rp ls 和 
rs 分 别 表示 左 行 、 右 行 了 波 和 S 波 。 时 域 递归 分 析 
方法 需要 将 微分 方程 转变 为 差分 方程 ,因此 需要 将 
计算 时 间 步 设置 到 一 定 小 的 程度 以 确保 计算 的 精 
度 。Li 等 .中 利用 时 域 递归 分 析 方 法 提出 了 确定 填 
充 节 理 法 向 和 切 向 刚度 的 方法 。 

以 上 理论 方法 各 有 优 缺 点 , 且 它 们 只 适用 于 分 
析 应 力 波 穿 过 单条 或 平行 的 无 限 延伸 节理 的 情况 。 
| 
法 才能 得 到 结果 。 但 是 理论 分 析 方 法 获得 的 结果 可 
UI E 
1.3 应 力 波 通 过 交叉 节理 

当 岩 体 中 存在 交叉 节理 组 时 ,应 力 波 的 传播 过 
程 变 得 非常 复杂 。Chai 等 ”采用 理论 方法 分 析 
了 平面 P 波 通过 两 条 交叉 节理 的 透射 规律 。 由 于 问 
题 的 复杂 性 ,该 理论 研究 忽略 了 应 力 波 在 节理 交叉 
] 力 学 学 报 
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点 处 的 散射 效应 。 对 于 含 多 组 交叉 节理 组 的 岩 体 ， 
数值 模拟 方法 是 研究 应 力 波 传播 的 有 效 方法 。 在 有 
限 元 法 和 有 限 差分 法 中 ,通常 用 特定 的 节理 单元 来 
表征 节理 !%] 。Song 等 “用 有 限 元 法 分 析 了 应 力 波 
作用 下 含 交叉 节理 组 的 岩 质 边 坡 的 变形 特征 。Dai 
等 .m1 用 有 限 差分 法 分 析 了 爆破 应 力 波 作用 下 含 多 
个 交叉 节理 组 的 岩 质 边 坡 的 位 移 响 应 。 然 而 这 样 的 
节理 处 理 方案 只 适用 于 节理 数量 少 且 小 变形 的 情 
况 。 相 对 而 言 ,离散 元 法 在 处 理应 力 波 通 过 复杂 分 
布 节理 问题 上 具有 明显 的 优势 。 在 离散 元 法 
(DEM) 中 , 岩 体 被 视 为 块 体 的 组 合 ,而 节理 作为 块 
体 之 间 的 界面 。Chen 等 ”用 UDEC( 二 维 离散 
元 软件 ) 模 拟 了 爆破 产生 的 应 力 波 在 包含 两 条 正 交 
节理 的 岩 体 中 的 传播 。 Zhu 等 ”用 UDEC 分 析 了 应 
力 滤 穿 过 多 个 交叉 节理 组 的 情况 ,系统 地 探讨 了 节 
理 的 力学 和 空间 分 布 包 括 节理 法 向 刚度 和 剪 切 刚 
麻 节 理 间距 .节理 组 之 间 的 交叉 角 、 入 射 波 方向 和 
AE TCU RE. ZUG. Zhu 等 ”用 UDEC 分 析 
了 爆破 诱发 的 应 力 波 通过 含 交叉 节理 的 岩 体 后 引起 
的 地 下 洞 室 围 岩 破坏 情况 。 其 他 的 离散 元 法 ,如 不 
连续 变形 分 析 方法 (DDA) “有 限 元 -离散 元 耦合 
PE (FDEM) 7 也 被 开发 用 于 研究 应 力 波 在 复杂 着 
体 思 传播 规律 。 值 得 指出 的 是 ,目前 大 部 分 的 数值 
分 办 限于 应 力 波 二 维 传播 问题 。 然 而 二 维 模型 难以 
全 面 反映 三 维 空间 中 节理 组 分 布 情况 。Deng 467 
用 BPEC 分 析 了 地 下 洞 室 爆 破产 生 的 应 力 波 通过 3 
组 种 理 组 的 传播 规律 。Haghnejad 等 ”用 3DEC 分 
术 生 节理 组 的 倾向 对 应 力 波 传播 规律 的 影响 。Li 
等 Gy 用 3DEC 分 析 了 爆破 应 力 波 作用 下 两 组 节理 
的 岩 体 中 洞 室 围 岩 的 变形 特征 。 上 述 数值 研究 的 关 
键 是 要 首先 获取 可 靠 的 节理 力学 参数 ,从 而 为 准确 
预测 应 力 波 在 复杂 岩 体 中 衰减 规律 提供 保障 。 


2 岩 体 动态 等 效 连续 介质 方法 


岩 体 动态 等 效 连续 介质 方法 是 分 析 应 力 波 在 兰 
体 中 传播 规律 的 另 一 有 效 方法 。 等 效 介 质 法 将 岩 体 
视 作 连续 介质 ,用 等 效 介质 的 有 效 模 量 来 描述 节理 
对 岩 体 力学 特征 的 影响 。 节 理 通常 以 平行 成 组 的 形 
SU ,导致 岩 体力 学 性 质 呈 现 各 向 异性 。 典 
型 的 含 一 组 平行 节理 的 宕 体 等 效 模 量 计算 模型 如 图 
6 所 示 。 通 过 5 个 弹性 常数 来 描述 该 岩 体 横 观 各 向 
同性 ,其 应 力 和 应 变 关系 为 
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K6 含 一 组 节理 的 岩 体 等 效 模 量 计算 模型 


Fig.6 Calculation model of equivalent modulus 


of rock mass containing a set of parallel joints 
基于 准 静 态 或 长 波长 假设 ,文献 [81- 82] 采 用 
平均 应 变法 推导 了 这 5 个 弹性 常数 与 完整 岩石 弹性 
参数 节理 间距 和 节理 度 之 间 的 关系 , 即 
Cy =p/(y +E,) 
Cy =p/(1 + £,) 


pee (18) 
DoD ey 
Cj, [1 *4(1-y)E,] C. 
C, = Cy, -2Cq 
其 中 
sacks) 
E, u/ (Sk,) (19) 


y (1-2v)/(2 -2v) 

SUP us 为 完整 岩石 的 剪 切 模 量 ;v 为 完整 岩石 的 泊 
松 比 ;S 为 节理 间距 ;到 A k, 分 别 为 节理 的 切 向 和 法 
I] ELI RE ; E, 和 万 , 分 别 表 示 节 理 对 岩 体 力学 性 质 的 
贡献 , 当 后 co 一 oo 时 ,有 cs 一 0 ,此 时 岩 体 的 弹性 常数 
等 于 完整 岩石 的 弹性 参数 。 

Daley 和 Thomsen ^! 给 出 了 了 波 、SV 波 和 SH 
波 在 横 观 各 向 同性 介质 中 的 相 速度 表达 式 为 


1 
pc (0) 7516s *C4*(C, - C4)sin'0 € D(0)] 


pé (0) = 六 [Ca + Cu + (Cu 7 Cs )sin'8 - DCA)] 


pcs (0) = Ca sin^8 + C,,cos^0 
(20) 
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其 中 :9 为 应 力 波 传播 方向 与 节理 法 向 之 间 的 夹 角 ; 
CP、C csv 和 cs 分别 为 P 波 SV 波 和 SH 波 波 速 。 
D(@) = {( C3; E Cay +4 ( C + Cm) sin 0 一 
2[ (C4 - C4) (Cy, + C4 -2C,4) ]sin^8 + 
[CCa 4045-20.) -4 (C4 Cay Jem B) ^ 
(21) 
由 于 上 述 等 效 介质 法 将 含 节理 的 岩 体 等 效 为 理 
想 弹 性 介质 , 且 无 法 考虑 应 力 波 在 节理 之 间 多 次 透 
反射 情况 ,因此 该 方法 难以 考虑 波 速 的 频率 相关 性 
和 透射 波 的 幅 值 衰减 规律 Li 等 ”指出 等 效 介 
质 法 用 等 效 模 量 隐 式 地 考虑 了 节理 的 作用 ,忽略 了 
节理 组 对 应 力 波 的 多 次 透 反 射 效应 。 
Pyrak-Nolte ^ BU Fi $t f P E LA] I 
力 波 在 岩 体 中 传播 时 出 现 的 频率 相关 的 波 速 变化 和 
赋值 误 减 效应 。 乎 弹性 模型 可 以 反应 固体 的 黏 弹性 
考 性 ,可 用 于 分 析 应 力 波 在 黏 弹性 介质 中 的 传播 规 
em. Ut B3 $5 sd PEU R15 Kelvin 模型 |、 
Maxell BUS VL Be J7 3 HU Bist pego m, Li agen 
BI" NANYE , BN Voiget KAIR Sc PES IE 
形成 的 三 单元 体 (图 7) 分 析 了 频 域内 平面 P 波 垂直 
完了 过 一 组 平行 线 ; 性 节理 的 衰减 规律 (图 8a) ,并 引入 
RESO BE ( FA 8c) 考 虑 了 应 力 波 在 节理 之 间 的 多 次 
ad. VAJ LA TA ERE U E E A ^E PI A 
% E A 


(3 +E Jo +n, 2 


tb E, WAT RERI E, 为 节理 所 贡献 
MISERERE n, WEER 
L : 


nE, OE Ee = 0 (22) 


1, 
图 7 Voiget 体 和 弹簧 元 件 串联 的 三 单元 体 


Fig.7  Three-element viscoelastic model composes 


of a spring in series with a Voigt model 
对 于 了 波 的 一 维 传播 问题 , 式 (22 ) 的 质点 振动 
速度 的 通 解 为 

v=A,exp(Bx)exp[i( wt — ax)] (23) 
其 中 :4 为 人 射 波 的 幅 值 ;w 为 圆 频率 ;B 为 应 力 波 
幅 值 衰减 系数 ;a 表示 单位 长 度 的 相位 偏 移 量 。 根 
据 位 移 不 连续 方法 , 简 谐 P X v, = A, exp (iot) 38 
单条 节理 后 的 透射 波 为 


Chinay n (EB HH TI 
dg CI INaA IVE TF RAT | 第 39 卷 
2k./z, 


w +(2k,/z,)° 


Vy = Ån X exp[ i( wt —xw/c, +t,)] 


(24) 

I -@/(2k,/z,)] (25) 

当 节 理 间距 为 5, 由 节理 和 岩石 组 成 代表 性 单 

位 体 的 长 度 为 5, 因 此 在 x=5 处 ,由 式 (23) 得 到 的 

质点 振动 速度 应 该 等 于 由 式 (24) 得 到 的 结果 ,由 此 
求解 得 到 式 (23) 中 的 未 知 参数 6 和 a 为 
2k,/z, | 


1 
n sn| for + (2k,/z,)° 
o 1 af o 
E TEM Era /z, ] 
MX (Q3) n] DUI, SCR SEH J JO ERU ER PA 
可 以 考虑 波 速 频 散 效应 和 幅 值 衰减 效应 ,但 是 忽略 
了 应 力 波 在 节理 之 间 的 多 次 透 反 射 效 应 。 为 了 解决 
该 问题 ,假设 每 个 节理 对 应 的 位 置 处 都 有 一 个 虚拟 
波源 ,其 作用 是 当 任 何 入射 透射 或 反射 波 到 达 节 理 
位 置 时 ,都 会 有 应 力 波 从 虚拟 波源 产生 。 虚 拟 波源 
产生 的 应 力 波 应 该 确保 整个 系统 能 量 守恒 , 即 人 射 
波 产生 的 能 量 等 于 通过 该 岩 体 的 透射 波 和 反射 波 产 
生 的 能 量 之 和 。 例 如 当 入 射 P 波 从 a nidi SIS b 
点 时 ,根据 式 (23) ,5b 点 的 质点 振 速 为 
v.(t,S) =A,exp(BS)exp[i(wt -a$S)] (27) 
adds e 表示 基于 等 效 黏 弹性 介质 模型 。 根 据 能 
守恒 定理 和 Kramer-Kronig K A, b 点 处 由 虚拟 波 
源 产生 的 反射 波 为 
v.(t,5) 2A, V1 - [exp(gS)]^ x 
exp[i(a@t + aS — 17/2) ] (28) 
此 时 ,点 处 的 入 射流 和 虚拟 波源 产生 的 反射 
波 沿 着 相反 方向 分 别传 播 至 a 和 c 点 位 置 ,然后 ,a 
All 点 的 虚拟 波源 产生 新 的 反射 波 ,并 周而复始 。 
因此 ,通过 等 效 黏 弹性 介质 的 透射 波 (图 8 P d A) 
为 不 同时 刻 到 达 的 多 个 透射 波 的 车 加 。 
Li 等 ' 引 证 明了 此 方法 获得 的 透射 波 与 基于 特 
征 线 方法 得 到 的 结果 一 致 。 之 后 , Ma S 用 相同 
的 方法 分 析 了 了 波 垂直 通过 多 条 平行 大 间距 节理 时 
的 速度 频 散 效应 和 幅 值 衰减 效应 。Fan 等 “提出 了 
考虑 节理 非 线 性 变形 特征 的 岩 体 等 效 昔 弹性 介质 模 
型 ,用 以 分 析 P 波 垂直 通过 等 间距 非 线 性 节理 的 透 
射 规律 。 值 得 注意 的 是 上 述 等 效 黏 弹性 介质 法 局 限 
于 解决 应 力 波 一 维 传播 规律 的 问题 。 当 岩 体 中 存在 
多 组 节理 时 , 岩 体 可 视 为 各 向 同性 介质 ,Li EO OR 


t, = tan 


(26) 
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用 Kelvin 模型 并 结合 岩 体 动态 品质 因子 ,得 到 各 向 同 
性 岩 体 的 等 效 连 续 介 质 模型 ,并 在 此 基础 上 分 析 应 力 
波 通过 多 层 岩 体 的 传播 问题 “””) ,该 方法 为 应 力 波 
传播 规律 分 析 方 法 的 实际 工程 应 用 提供 了 基础 。 


—— | 


1 
1 
站 L ee en apes i PE ae tas ma EUM TP] (omms IER 
oe x j z 
S S $ 
m 
> e 
(a) 含 等 间距 节理 的 岩 体 
(ui pehrseseed 
CO! | 
e | | ' 
I 1 1 1 
e | | | 
iJ i i i 
mi a earl peers as | 1 
Ts E 
N 
N (b) 代表 体 单元 (c) 图 (b) 的 等 效 单元 体 
e VWS VWS VWS 
QU | | ep i eas ee als aa ay) i 
mu Y | | | | | 
2 | | es: l 1 EN 
» | [Eee | ! 
a: l l I I i 
[Ph T eee ST (ae T EE a Te E 
c E 5 E d 
E ans S S S 
Q 


(d) 含 虚拟 波源 的 等 效 介质 模型 
图 8 基于 虚拟 波源 的 岩 体 等 效 黏 弹性 介质 模型 方法 
Fig.8 Equivalent viscoelastic medium model 


. n ti » 
including “virtual wave source 


3 双 斥 度 非 连续 耦合 方法 

随 着 极端 环境 下 岩 体 工程 的 发 展 , 高 低温 岩 体 
力学 特性 日 益 受 到 重视 。 在 温度 的 影响 下 , 岩 体 内 
出 现 不 可 忽略 的 细 观 裂 际 ,如 图 9 所 示 。 温 度 导致 
的 岩 体 内 细 观 裂隙 的 产生 会 进一步 影响 应 力 波 的 传 
播 。 含 细 观 裂 阶 的 岩 体 通常 被 等 效 为 完全 弹性 的 介 
质 , 采 用 等 效 弹性 模型 研究 应 力 波 的 传播 特性 。 男 
外 ,实验 研究 表明 应 力 波 在 含 细 观 裂隙 的 岩 体内 传 
播 时 会 发 生 速度 频 散 效应 和 幅 值 衰减 效应 ””。 
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Johnston ^ 和 刘建华 等 ”将 含 细 观 裂隙 的 岩 体 内 波 
的 衰减 归结 为 由 摩擦 ,流体 流动 黏 性 松弛 和 散射 造 
成 的 。 因 此 ,等 效 黏 弹性 介质 被 提出 用 于 研究 含 细 
观 裂 隙 的 岩 体 的 动态 力学 行为 ”| Ren 等 ' 基于 
数值 模拟 的 方法 进一步 验证 了 含 细 观 裂隙 的 宕 体 的 
等 效 黏 弹性 特性 。Fan 4E?" 采用 Maxwell 体 和 弹 得 
元 件 并 联 的 三 单元 体 来 描述 含 细 观 裂 院 的 岩 体 的 黏 
弹性 行为 。Niu 25 7" 通过 数值 和 试验 相 结合 的 方 
法 Wine T E MARR A AY BE A BK, 


pe: ede 
ee 


A9 岩 体 内 的 细 观 裂隙 
Fig.9 Micro-defects in rock mass 
基于 含 细 观 裂隙 的 岩 体 的 等 效 忒 弹性 模型 ,三 
特征 线 方法 被 提出 来 研究 应 力 波 在 竹 弹 性 介质 中 的 
传播 (图 10) 。 根 据 应 力 波 传播 的 运动 方程 为 


p® - 88.9 (29) 
其 中 p 为 岩 体 密度 。 
应 力 波 传播 的 连续 方程 为 
BE e (30) 
ðt Ox 


255 ARE Ap Jot Bie A BY AS 9475 Be, n] VAR — 
族 特征 线 和 对 应 的 相 容 关 系 。 从 而 建立 计算 含 细 观 
裂隙 的 岩 体 内 应 力 波 传播 的 三 特征 线 方法 ,如 图 10 
所 示 。 根 据 三 族 特 征 线 和 特征 线 对 应 的 相 容 关 系 可 
以 求 得 含 细 观 裂 际 的 岩 体 内 应 力 波 传播 的 质点 应 
力 \ 质 点 速度 和 质点 应 变 。 


图 10 三 特征 线 方法 


Fig. 10 Three characteristic lines method 
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细 观 缺陷 节理 


图 11 含 双 尺度 结构 面 岩 体 的 处 理 方案 
Fig.11 Scheme of modeling of rock mass 
with double-scale discontinuities 

C 2 力 波 的 传 
播 耻 存在 显著 的 影响 。 不 同 尺度 的 不 连续 结构 面 影 
S ATE TEL: a. cU LLG 
VIE ce PR D ERE FR SS AZ ED REA 
细 观 裂 阶 使 得 应 力 波 产生 弥散 和 衰减 。 因 此 ,建立 
统一 模型 描述 含 细 观 裂 阶 和 宏观 节理 的 岩 体 ,并 在 
Bes HE ae ec pec pe ellie 
PSR ES FEE FE TAMER. PAG, Fam EUO 将 等 效 介 
B 1 ui us 
PA MALWARE, C CRURA E T — 
Als TI We LE DUN BE TRE S rp D FE AEE 
TOT Gr AA RLR H A FG BON BH FEE 
SD US ER OLN BAAS YE B WARE. EL 11 表 
ARP WR EAH THA PR HU AB TRU R, ER 
用 相同 的 处 理 方法 在 频 域内 计算 了 平面 了 波 以 任意 
角度 倾斜 通过 含 双 尺 度 结构 面 的 岩 体 的 透射 系数 。 


ARCU PE BUR 
CHENA) 


宏观 节理 


菱形 单元 
P GDP J) 


PAREEN 


分 离 式 三 特征 线 方法 
Fig. 12 Split three characteristic lines method 


针对 含 双 尺 度 结构 面 岩 体 的 均匀 不 连续 建 模 思 


图 12 
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路 ,Fan 等 "将 三 特征 线 与 位 移 不 连续 相 结合 , 提 
出 了 分 离 式 三 特征 线 方法 ,在 时 域内 分 析 了 应 力 波 
通过 含 双 尺度 结构 面 的 岩 体 的 透射 系数 和 有 效 速 
度 。 分 离 式 三 特征 线 方法 如 图 12 所 示 。 与 传统 的 
三 特征 线 方法 不 同 ,分 离 式 三 特征 线 方法 中 除了 计 
算 边 界 点 和 岩 体 内 部 点 处 应 力 波 传播 的 三 角形 单元 
和 菱形 单元 ,还 引入 了 计算 应 力 波 在 宏观 节理 处 传 
播 的 分 离 式 次 形 单元 。 分 离 式 三 特征 线 方法 中 的 三 
种 基本 单元 ,如 图 12 所 示 。 以 上 研究 均 基 于 宏观 节 
理 为 线性 变形 的 情况 ,适用 于 宏观 节理 位 于 应 力 波 
远 场 时 。 然 而 , 当 宏 观 节理 位 于 应 力 波 近 场 或 者 合 
有 填充 物质 时 , 需 考虑 宏观 节理 的 非 线 性 变形 。 针 
对 应 力 波 穿 越 非 线性 宏观 节理 的 情况 ,Fan 等 "中 
通过 引入 分 段 线性 位 移 不 连续 方法 ,改进 了 王 礼 
立 5 史 提出 的 三 特征 线 时 域 分 析 方法 ,分 析 了 非 线 性 
PE E A E D a 
的 影响 。 随 后 ,通过 将 多 脉冲 作用 下 宏观 节理 的 力 
es 
冲 应 力 波 在 含 双 尺度 结构 面 的 岩 体 内 的 传播 特 
PE) 。 另 外 ,Huang E°! 通过 改进 矩阵 传递 法 在 
频 域内 计算 了 应 力 波 在 倾斜 穿越 含 双 尺 度 结构 面 的 
岩 体 时 的 透射 .反射 系数 。 
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岩 体 中 广泛 存在 的 节理 对 应 力 波 的 幅 值 衰减 和 
速度 频 散 具有 明显 的 影响 。 近 年 来 ,学 者 们 提出 了 
不 同 的 理论 分 析 方 法 用 于 分 析 应 力 波 在 岩 体 中 的 传 
播 规律 研究 。 本 研究 对 现 有 方法 进行 了 详细 地 综 
述 ,并 探讨 了 它们 的 适用 范围 。 

在 非 连续 介质 方法 中 , 基于 位 移 不 连续 方法 的 
特征 线 法 、 时 域 递 归 分 析 法 .和 矩阵 传递 法 已 被 广泛 地 
用 于 分 析 应 力 波 通过 单条 和 一 组 节理 的 问题 。 相 对 
于 频 域 分 析 方 法 ,时 域 分 析 方 法 由 于 可 直接 分 析 任 
意 波形 的 人 射 波 的 传播 问题 ,从 而 简化 求解 过 程 。 
然而 ,位 移 不 连续 方法 不 适用 于 应 力 波 通过 交叉 节 
理 组 的 复杂 情况 。 已 有 的 数值 方法 ,如 有 限 元 法 、 有 
限 差分 法 离散 元 法 以 及 耦合 方法 可 以 分 析 应 力 波 
通过 节理 复杂 空间 分 布 的 情况 ,但 是 它们 面临 的 两 
个 关键 问题 ,首先 是 节理 的 动态 力学 模型 准确 描述 ， 
其 次 ,数值 算法 有 竺 改进 从 而 提高 计算 性 能 。 

对 于 动态 等 效 连续 介质 法 ,目前 提出 的 考虑 虚 
拟 波 源 的 岩 体 条 弹 性 动态 等 效 连 续 介质 方法 能 够 解 
力学 学 报 
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释 应 力 波 在 含 结构 面 的 岩 体 中 传播 时 发 生 的 频率 相 
关 的 幅 值 衰减 和 速度 降低 的 现象 ,但 是 该 方法 局 限 
于 应 力 波 一 维 传播 问题 ,对 于 高 维度 的 问题 ,目前 在 
一 些 特定 情况 下 提出 了 基于 岩 体 动态 品质 因子 的 黏 
弹性 等 效 介质 方法 ,该 方法 还 有 待 开 展 进一步 的 
研究 。 

连续 和 非 连 续 耦 合 方法 适用 于 分 析 岩 体 中 含有 
宏观 节理 和 细 观 千石 裂 际 ,该 方法 目前 侧重 于 应 力 
波 一 维 传播 规律 研究 ,如果 和 基于 岩 体 动态 品质 因 
子 的 黏 弹性 等 效 介质 方法 相 结合 ,可 开展 应 力 波 二 
维 传 播 问 题 的 研究 。 
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